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ABSTRACT 

Black hole masses in active galactic nuclei (AGN) are difficult to measure using conventional dy- 
namical methods, but can be determined using the technique of reverberation mapping. However, it 
is important to verify that the results of these different methods are equivalent. This can be done in- 
directly, using scaling relations between the black hole and the host galaxy spheroid. For this purpose, 
we have obtained new measurements of the bulge stellar velocity dispersion, cr*, in Seyfert 1 galaxies. 
These are used in conjunction with the M, - cr* relation to validate nuclear black hole masses, M,, in 
active galaxies determined through reverberation mapping. We find that Seyfert galaxies follow the 
same M, - cr* relation as non-active galaxies, indicating that reverberation mapping measurements 
of M, are consistent with those obtained using other methods. We also reconsider the relationship 
between bulge absolute magnitude, M^ui, and black hole mass. We find that Seyfert galaxies are offset 
from non-active galaxies, but that the deviation can be entirely understood as a difference in bulge 
luminosity, not black hole mass; Seyfert hosts are brighter than normal galaxies for a given value of 
their velocity dispersion, perhaps as a result of younger stellar populations. 

Subject headings: galaxies: nuclei — galaxies: Seyfert — quasars:general — galaxies: kinematics and 
dynamics 



1. INTRODUCTION 

An observational problem in research on active galac- 
tic nuclei (AGN) is determining the masses of the cen- 
tral black holes. Application of stellar and gas dy- 
namical tedungu^s-conm galax- 
ies fe.g. lBower eTaLHTool ICebhardt et al.ll2003j) . is dif- 
ficult or impossible in type 1 AGN due to the bright 
nucleus. A potential solution is the technique known 
as reverberation mapp ing (Blandford & Mckec 1982; 
iNetzer fc Petersoiilll997j) . which has been used to esti- 
mate black hole masses for a bout 30 Seyfert gal axies and 
quasars ijWandel et al.iri999HKaspi et alJl200(]|) . Rever- 

^ Visiting Astronomer, Kitt Peak National Observatory, Na- 
tional Optical Astronomy Observatory, which is operated by 
the Association of Universities for Research in Astronomy, Inc. 
(AURA) under cooperative agreement with the National Science 
Foundation. 

^ Kitt Peak National Observatory and National Optical Astron- 
omy Observatory are operated by the Association of Universities 
for Research in Astronomy, Inc. (AURA) under cooperative agree- 
ment with the National Science Foundation. 



beration mapping allows the radius of the broad line re- 
gion, Rblr, to be determined by attributing the time 
lag between variations in the continuum and line emis- 
sion to light travel time effects. Assuming that the ve- 
locity width, VF, of the permitted emission lines is due to 
Keplerian motion in the black hole potential, the mass 
can estimated by M, ~ W^Rblb/G, where G is the 
gravitational constant. Thus, this method may provide 
an opportunity to study the role of the nuclear gravi- 
tational potential in large numbers of AGN. But it has 
not yet been demonstrated that reverberation mapping 
yields Af, values that are consistent with those obtained 
from analysis of the stellar and gas dynamics. Given the 
difficulty of a direct approach, indirect methods must be 
used to compare these different techniques. Our purpose 
in this study is to compare reverberation mapping and 
stellar dynamical measures of black hole mass using scal- 
ing relations between the black hole and the host galaxy 
bulge. 

The correlation between black hole mass and 
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bulge absolute maKnitude, Mh,d, ijKormendv I 119931: 
iKormendv fc Richstone |I199,'tD . can be used to compare 
results of these different approaches. However, compar- 
isons of the M,-Mbui relation in active and normal galax- 
ies have produced conflicting results. Two early studies 
||Ho,1999.: ..Wandel 1999) reported that the reverberation 
mapping masses for AGN were smaller than the stellar 
dynamical masses in normal galaxies of the same bulge 
luminosity. This called into question the validity of the 
reverberation mapping masses, suggesting a systematic 
trend to unde restimate M, b y a factor of ~ 5 (,Ho„,199i^ . 
Alternatively, iWandell l)1999D suggested that the masses 
were accurate and that there were real differences in the 
black hole masses of active and non-acti v e gala xies of 
equal bulge mass. More recently, 'Wandcl ( 2002'), revis- 
ing his early assessment, and McLurc & Dunloo (200l|), 
using bulge magnitudes determined in most of their sam- 
ple from Hubble Space Telescope (HST) images and as- 
suming disk-like kinematics for the BLR, find that active 
galaxies follow the same M, - M^ui relation as normal 
galaxies. 

It is difficult to make a fair assessment of these differ- 
ing results since some of these studies have used bulge 
magnitudes derived from the relation between bulge-to- 
disk ratio as a function of morphological type formulated 
by Simicn & dc Vaucoulcurs (1986), while others have 
applied bulge-disk decomposition techniques to ground- 
based or HST images. Furthermore, bulges of Seyfert 
galaxies tend to be small — their typical effective radii 
in all but the nearest objects are a few arcseconds or 
less. This, combined with the effects of the bright nuclear 
source in type 1 objects, which are of course the candi- 
dates for reverberation mapping, makes ground based 
estimates of M^ui quite tricky. 

The relationship between black hole mass and bulge 
stellar veloc ity dispe rsion^ cr*, (|Ferrarese & MerritO 
l200rllc"ebhardt et a l. 2000a), provides an alternative and 
independent way to compare stellar dynamical and re- 
verberation mapping results. In this relation, which is 
stronger than the one using Mfe„j , cr* is measured in aper- 
tures extendin g to distances of a kiloparsec or so from the 
nucleus ( Gebhardt et alJl2000al use apertures extending 
to the bulge effective radius, r^, and denote their mea- 
surement CTe)- The goal is to standardize the region over 
which the dispersion is determined and more importantly 
to minimize the influence of the nuclear black hole on the 
stellar kinematics. 

The small number of published cr* measurements for 
type 1 AGN, has until re cently prevented a stati stically 
significant investigation. IGebhardt et al.l l|2000bj) found 
no deviation for 7 Seyfert Is from th e M, — t Tf. re lation. 
The same conclusion was reached bv iNelsonl f2000), us- 
ing FWHM[oiii] in place of ct*, on the assumption that 
the forbidden line kinematics in Seyfert galaxies are pre- 
dominantly the result of virial motion in t he bul ge grav- 
itational potent i al (e. g. |Nelson & Whittll [T996j) . Also, 
IFerrarese et al.l 11200 11 obtained tr, measurements for 6 
Seyfert 1 galaxies and find these values to be consistent 
with expectations from the M, — ct, relation. 

In this paper we present new measurements of cr, in 14 
Seyfert 1 galaxies to check the M, values obtained from 
reverberation mapping for consistency with the M, — 
relation established by stellar dynamical techniques. Ob- 
servations of the Ca II triplet lines and velocity disper- 
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Velocity Dispersions in Seyfert 1 Galaxies 



Name 


Aperture 
\") 


(T»±e 


C2±e 

I iS.lil / o j 


log m; 

i Mr-, \ 


NGC 3227 


4.7 


136±4 


1187±5 


r.OD_Q 21 


NGC 3516 


4.3 


181±5 


2657±6 


7 r,q + 0.08 

'•^•^-0.09 


NGC 4051 


5.6 


89±3 


743±6 


6 15+0-32 
"•13-0.45 


NGC 4151 


6.5 


97±3 


1020±6 


7-08lS:i 


NGC 4593 


6.5 


135±6 


2602±7 


6-82l°:25 


NGC 5548 


2.2 


201±12 


5100±15 


7 Q7+0.07 
'•="-0.07 


NGC 7469 


4.7 


131±5 


4947±7 


6.88l°:i» 


Mrk 79 


4.1 


130±12 


6540±8 




Mrk 110 


4.7 


91±25 


10509±13 


6.89«:» 


Mrk 279 
Mrk 335 
Mrk 509 


4.7 


197±12 


9062±12 


6.58«:1^ 
7 96+0-05 


Mrk 590 


4.7 


189±6 


7904±7 


7 14+0.10 


Mrk 817 


4.7 


120±15 


9409±11 


7 c:c:+0.11 


Akn 120 


3.0 


221±17 


9575±15 


*-^'-0.12 


3C 390.3 


3.6 


273±16 


16640±15 


o c-7+0.12 
».tlC_0 21 



* The black hole masses are th ose using the r.m.s. BLR widths. 
Sources are IQnken ct alj l2()(1^ for NGC 3227, NGC 3516 and 
NGC 4593, iHo 11999) for NGC 7469 and IKaspi et all l200Sh for 

the rest. 



sion measurements using the Fourier Correlation Quo- 
tient (FCQ) are presented in section|21 Bulge magnitude 
estimates for 11 of these obtained from Hubble Space 
Telescope archival data are presented in sectional Notes 
on individual sample members are presented in section 
01 In section Owe discuss the results in the context of 
the M, — (7* relation and in section we revisit the Af, 
- Mb„; relation to try to determine the origin of the dis- 
crepancies noted above. Our results are summarized in 
section 

2. SPECTROSCOPIC DATA 

2.1. Observations and Reductions 

Longslit spectroscopy of 16 Seyfert 1 galaxies and sev- 
eral stellar template stars was obtained on two observ- 
ing runs (U.T. dates April 10-14, 2001 and November 
1-4, 2001) using the 4-meter telescope at Kitt Peak Na- 
tional Observatory. Our sample is composed of Seyfert 1 
galaxies included in long term emission line and contin- 
uum monitoring campaigns for which reverb eration map- 
ping estimates of M , have be en pu blished ijKaspi et alJ 
l200nt IWandel et all 1999; Hol lT999l) . The RC spectro- 
graph was used with the LBIA CCD which has high 
quantum efhciency in the far-red ( ~ 90% at 8500 A), 
excellent cosmetics and no fringing l)Dev e t al."200li) — 
an ideal detector for observing the Ca II triplet lines 
(A8498.0, A8542.1 and A8662.1 A). These properties al- 
lowed us to use a high dispersion grating (BL 380 giving 
0.46 A/pixel in the raw data) and a narrow slit (1" ) to 
achieve high spectral resolution (~ l.lA or 40 km 
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Fig. 1. — Reduced spectra are presented with telluric absorption correction and relative flux calibration applied. The position of the 
Ca II triplet lines are indicated by the bar with three ticks below each spectrum. The tick marks on the left hand vertical axis indicate the 
zero points for each spectrum, giving a sense of the broad range in Ca II triplet equivalent width. 
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Fig. 2. — Wc show the finished spectra which were input to FCQ and also the broadened stellar template. The spectra have been scaled 
so that the middle line in the Ca II triplet has approximately the same depth in each galaxy. 
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FWHM). We used an order-blocking filter (OG570) to 
prevent contamination by second-order blue light, which 
will of course be strong in Seyfert 1 spectra. Exposure 
times ranged between 3600 seconds for nearby bright 
galaxies (e.g. NGC 4051, NGC 3516) and 10,800 sec 
for fainter objects (e.g. Mrk 110). Standard spectro- 
scopic reductions — bias subtraction, flat field correc- 
tion, geometric distortion and wavelength calibration — 
were carried out using IRAF. Conditions were generally 
not photometric so only a relative flux calibration was 
applied. 

The strong telluric emission and absorption features 
in the far-red present a difficult reduction problem, par- 
ticularly for the higher redshift galaxies in our sample. 
For targets with z > 0.016, where the Ca II triplet 
and night sky emission lines begin to overlap (except 
3C 390.3, see section QJ, sky-subtraction was accom- 
plished usin g the "nod- and-shufile" observation mode, 
developed bv lGlazebrook & Bland-Hawthorn (2001) and 
recently implemented at Kitt Peak. Briefiy, by repeat- 
edly re-positioning the telescope (nodding) and shifting 
the charge on the detector (shuffling), both the target 
and the sky are observed at the same spot along the 
slit and in the same pixels on the detector, permitting 
excellent sky subtraction. The drawbacks for the nod- 
and-shuffle technique are a factor of ^2 increase in the 
sky noise and extra overhead time for the telescope nods. 
Although there were difficulties in the implementation of 
this mode at the time, we were quite pleased with the re- 
sults and we recommend use of the perfected observing 
mode where sky subtraction is critical and signal-to-noise 
is not of primary concern. 

Atmospheric absorption bands begin to be a problem 
at wavelengths above 8900 A, again affecting the higher 
redshift members of the sample. Even if the Ca II triplet 
lines are not directly affected, correction for these can 
be important for achieving a good fit to the continuum 
on the red side. We therefore obtained spectra of white 
dwarf spectroscopic standard stars, matched in airmass 
and observation time to our Seyfert exposures, to use 
as atmospheric templates. Experimentation showed that 
scaling the equivalent width of the atmospheric bands, 
by adding a constant value to the template and renor- 
malizing, improved the results in some cases, helping to 
recover good line profiles for the reddest Ca II triplet 
line. In the case of 3C 390.3, the atmospheric absorption 
correction was quite critical since all of the Ca II triplet 
lines lie in a region of strong atmospheric absorption. In 
part due to the broad Ca II triplet features in this object, 
we were able remove the telluric lines successfully. 

Reduced, sky-subtracted nuclear spectra corrected for 
telluric absorption are shown in figure Q The spec- 
tra have been normalized to the mean counts in the 
spectral region near the Ca II triplet lines, whose po- 
sitions are indicated for each object. It is worth em- 
phasizing that a variety of spectral characteristics are 
present in these galaxies. The Ca II triplet lines show 
a range of equivalent widths presenting a challenge for 
analysis using the FCQ technique in some cases. Also 
a number of nuclear emission features, including broad 
01 A8446, [Fe II] A8616, Paschen and Ca II triplet it- 
self in emission, appear in our spectra and must be re- 
moved before analysis with FCQ. In fact for two galax- 



ies, Mrk 335 and Mrk 509, no velocity information 
could be obtained (see section^. Given the associa- 
tion of AGN and starbursts we must also keep in mind 
that the strength of the Ca II triplet lines may be en- 
hanced by young stell ar popula tions ric h in late-type 
supergiant stars (e.g. tTerlevich. Diaz fc TterlevichllT990l: 
INelson fc Whitid3l200ir 

2.2. Removal of Nuclear Emission 

Even in the far-red, the strength of the nuclear emis- 
sion from type 1 AGN relative to the contribution from 
the host galaxy makes obtaining high signal-to-noise 
spectra of the stellar absorption features a difficult task 
— effectively, the nuclear continuum is a source of noise. 
Also broad and narrow line nuclear emission, even when 
not affecting the Ca II triplet lines must be carefully re- 
moved before the velocity dispersion measurements can 
be made. In some cases, where the stellar absorption 
lines are strong and the nuclear emission is weak, it 
was possible to include the nucleus of the galaxy in our 
longslit extractions without further processing. For other 
galaxies, however, the central rows of the longslit spectra 
had to be corrected for contaminating emission features, 
or when the absorption features were entirely lost in the 
noise, excluded altogether. 01 A8446 and [Fe II] A8616 
could often be removed by subtraction of a single or dou- 
ble Gaussian fit. In other galaxies the blending of emis- 
sion lines and continuum was more complex, producing 
an irregular continuum shape near the Ca II triplet lines. 
In several cases this required subtraction of a model of 
the unresolved nuclear spectrum prior to the dispersion 
measurements . 

To produce the model, a narrow spectrum, usually two 
or three pixels wide centered on the nucleus was ex- 
tracted. Next, the weak Ca II triplet absorption fea- 
tures were removed from the model by subtracting the 
spectrum of a stellar template, which was broadened and 
scaled until the absorption lines vanished. Then, using 
spectra of nearby ninth or tenth magnitude SAO stars 
obtained prior to each galaxy exposure, we produced a 
seeing profile along the slit. This was combined with the 
model nuclear spectrum and subtracted from the origi- 
nal longslit data to correct for the contamination by the 
nuclear emission. The most important effect of this pro- 
cedure is that it removes the irregular continuum shape 
from spectra immediately adjacent to the nucleus. This 
was particularly helpful in NGC 4051, NGC 4151 and 
NGC 4593, where the continuum on the blue side of the 
spectrum is distinctly uneven (compare figures d and EJ . 
This allowed us to extract spectra for the velocity dis- 
persion analysis much closer to the nucleus than would 
have been possible otherwise, resulting in higher signal- 
to-noise and ensuring that as much of the bulge light was 
included as possible. 

2.3. Velocity Dispersion Measurements 

Velocity dispersions were obtained by analysis of the 
bulge line-of-sight velocity distribution (LOSVD) using 
the Fourier-Correlation Quotient method (Be n de r 19 9(|). 
For each galaxy a spectral region, including only the 
Ca II triplet lines and a small section of continuum on 
either side, was used in the analysis (rest wavelength 
8460 A - 8716 A). In most cases, the width of the spec- 
trum along the slit was chosen to cover at least 1 kpc 
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Fig. 3. — Our fTj, measurements are compared with those of 
IKerrarese eTall 120011) . The solid line is X = y. 



across each galaxy. For a few nearby galaxies we found 
that extending the slit all the way out to Ve (see table |3l 
actually degraded the signal-to-noise while yielding ef- 
fectively the same dispersion. In these cases we opted 
for the shorter slit length. For the fainter targets an 
optimal extraction maximizing signal-to-noise was used. 
Although this does not generally correspond to r^. as used 
by Gcbhardt ct al. ( 2000^ in defining the A/, — rela- 
tion, given the difhculties posed by the nuclear emission, 
it was a reasonable and unavoidable compromise. Nev- 
ertheless, our spectra include bulge starlight far enough 
from the nucleus to be beyond the influence of the black 
hole, while excluding emission from the kinematically 
colder galaxy disk. 

Gauss-Hermite polynomials were fitted to the LOSVD 
to determine the bulge systemic velocity, K, and the 
bulge stellar velocity dispersion, cr,. The uncertainties 
in Vif and cr* were estimated by running the FCQ algo- 
rithm on Monte Carlo simulations of the data, created 
using a stellar template convolved with the LOSVD and 
adding in noise to match that of the observed spectrum. 
The results for the sample are given in table ^ includ- 
ing M, values. Figure |21 shows the data used in the FCQ 
analysis shifted to the rest frame along with the template 
broadened by the derived LOSVD. 

Velocity dis persions for six of our s ample galaxies were 
published bv iFerrarese et aTl l|2001j) . A comparison of 
our results with theirs is shown in figure O We find 
that our values agree with theirs quite closely, with a 
marginal tendency for our dispersions to be larger at the 
low a^: end. This may be a consequence of our narrower 
slit (they used 2" ) which excludes more disk light — a 
result that may be important for systems with smaller 
bulges. 

3. BULGE MAGNITUDES 



To investigate the reported deviations of AGN from the 
M, - Mbui relation discussed in section ^ we have used 
archival HST images to measure bulge magnitudes. Ten 
of the Seyfert 1 galaxies in our sample were observed 
in th e program of WFPC-2 snapshots (iMalkan et al. I 
[19981) using the F606W filter and another, NGC 4151, 
was observed with the F547M filter. The main diffi- 
culty in determining the bulge magnitudes from these 
images was removing the point source nucleus, which is 
unfortunately saturated in the F606W images. To ac- 
complish this, mo del PSFs, cr eated using the Tiny Tim 
software package l|Krist lll99fT() . were shifted, scaled and 
subtracted from the images. The scale factor which pro- 
duced a minimum r.m.s. deviation in the residual image 
in an annulus incorporating the innermost unsaturated 
pixels, was chosen as an initial guess. Since this tends 
to overestimate the nuclear luminosity, further trial-and- 
error scaling of the PSF around this value allowed us 
to fine tune the results; the final best subtraction was 
judged by examining the removal of peculiar features 
of the outer PSF (rays and ring like structures) in the 
residual images. Typically we found that a range of 
scale factors produced acceptable results. The result- 
ing uncertainty in the bulge magnitudes, which given the 
subjective nature of the PSF scaling we uniformly take 
to be ±0™2, was folded into our error estimates. Be- 
cause the galaxy bulges are reasonably large in the im- 
ages (re ~ 50 — 100 pixels in most cases) our results do 
not depend critically on PSF subtraction. The purpose 
is mainly to remove the extended wings of the unresolved 
nuclear emission and recover as much of the inner bulge 
light as possible. For most objects, we were able to ex- 
tend our bulge profile analysis down to about half an 
arcsecond. 

Other complications included removing the cosmic ray 
hits from the images and accounting for emission lines 
from the Narrow Line Region (NLR). About 90% of the 
cosmic rays were removed automatically using a script 
written in IRAF. Of the remaining cosmic rays those 
close to the nucleus were removed by hand. The broad 
bandpass of the F606W filter includes the [0111] A5007, 
[Nil] A6548, A6584, and Ha emission lines which arise 
in the NLR and will not be removed with the point-like 
nucleus. Using the STSDAS task CALCPHOT, and the 
model Seyfert 2 galaxy spectrum provided by the SYN- 
PHOT package, we were able to estimate the NLR flux by 
scaling the model s pectru m to match the [0111] flux tab- 
ulated in lWhittlel l|1992ar . The [0111] fluxes are almost 
all rated quality "a" by Whittle indicating that they were 
obtained through large apertures and that there is good 
agreement when more than one measurement has been 
published. Thus we are confident that we have accounted 
for the bulk of the NLR emission. This provided a cor- 
rection to our bulge magnitude estimates and is listed in 
Table □ 

Ideally, two-dimens i onal d irect fitting techniques, such 
as that of lPeng et alJ (|2002'), and a less restricted model, 
e.g. a Sersic function, could and should be used to in- 
vestigate the bulge morphologies of these galaxies. We 
note in the descriptions of the individual galaxies below 
(section 0J that the bulge profiles of Mrk 590 and Mrk 
817 are better fit by exponentials than r^/^ laws. Also 
some display inner spirals, dust lanes, lenses and bars 
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Flc. 1. Luminosity profiles and bulge-disk decomposition are shown for Seyfert 1 galaxies from HST images. The emission from the 
nucleus has already been subtracted from the profiles. Note the exponential bulges in Mrk 590 and Mrk 817. A horizontal line shows 
regions included in the fitting, (if not present all of the data was used) gaps are usually associated with bars. 
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that further comphcate the analysis. However, account- 
ing for these difficulties and the uncertainties associated 
with NLR emission and saturated PSFs is outside the 
scope of our program. Thus, in this paper, we adopt ei- 
ther the r^/'* law or exponential fits to the azimuthally 
averaged luminosity profiles and use these to estimate 
bulge magnitude. 

We used the IRAF task ELLIPSE to produce lumi- 
nosity profiles from the image isophotes which were then 
fitted with the sum of an exponential or r^/^-law bulge 
and an exponential disk. In carrying out the fits it is 
important to avoid getting trapped in spurious minima 
of x^- The procedure we adopted to ensure meaningful 
results, was to fit the disk region first. Then holding the 
disk parameters fixed we included a fit to the bulge. As 
a final step, we fitted both bulge and disk simultaneously 
allowing all parameters to vary, using the parameter val- 
ues from the previous steps as the initial guess. The 
luminosity profiles and fits are plotted in figure 01 Note 
the presence of bars or inner ovals in a few galaxies (e.g. 
NGC 4593, Mrk 79, NGC 7469), seen as a run of excess 
light over the fit, between the bulge and disk. We gen- 
erally excluded these regions from the analysis. Bulge 
magnitudes were then estimated from the scale length 
and surface brightness parameters of the fits, assum ing 
Ho = 75 km s'^ Mpc'^ (as used bv lKasni et aUEmd to 
calculate black hole mass). Foreground extinction from 
the Milky Way has been removed, but no correction for 
dust within the host galaxy has been applied. The re- 
sults are given in table[21and notes on the characteristics 
of individual galaxies are given in section ^ 

We estimated the error in the bulge magnitudes using 
Monte Carlo simulations. Random noise with the same 
r.m.s. deviation as the fit to the luminosity profiles was 
added to the data and the profile was refitted. The re- 
sulting standard deviations in the fit parameters after 
30 trials was then propagated through to estimate the 
error in bulge magnitudes from the fitting. These were 
then added in quadrature to the uncertainty in the nu- 
clear magnitudes and adopted as the uncertainties listed 
in table 121 We also include an uncertainty of 15% in the 
distances of each galaxy to compute the final errors in 
the absolute magnitudes. 

Bulge magnitudes for se veral of the same Seyfer t 
galaxies were tabulated bv iMcLure fc Dunlop I l|200H) . 
some measured using the same HST images, while 
others were taken from the phot ographic study of 
iBaggett. Baggett fc AndersonI ||1998H. In fi gure 151 we 
plot our vahies against the lMcLure' fc Dimlop I l|20(m re- 
sults, where we have corrected their values to our choice 
of Hq and to the F606W bandpass (since they quote Mt,ui 
in the R bandpass). The solid line shows X — Y, and 
suggests no systematic trends, although the large scat- 
ter (r.m.s. =1™3, compared with the combined error 
estimate of 1™0) does indicate some disagreement in a 
number of galaxies. Although McLure & Dmilop (200ll) 
were careful to avoid saturating the nucleus in their own 
observations of quasar hosts, they do n ot discuss how 
they handled the saturated nuclei in the ijMalkan et al. I 
Il998j) snapshots. Also, in Mrk 817 we have a rather large 
correction for NL R emission, which to our kn owledge was 
not considered bv lMcLure &: Dunlop 1 1)200 Ij) . In any case 
the large scatter, while disappointing, is not far outside 



TABLE 2 



Bulge Magnitudes in Seyfert 1 Galaxies 



Name 


or Un 


Qj. 


■m-bui 








(mag.) 


'{") 


(mag.) 


(mag.) 


(mag.) 




{F606W) 




(F606W) (F606W) 


(F606W) 


NGC 3227 


18.31 


2.63 


13.24 


0.09 


-18.19 




±0.16 


±0.26 


±0.36 




±0.49 


NGC 3516 


19.43 


7.41 


11.91 


0.02 


-21.19 




±0.12 


±0.55 


±0.29 




±0.44 


NGC 4051 


18.58 


3.10 


13.08 


0.05 


-18.53 




±0.12 


±0.23 


±0.30 




±0.44 


NGC 4151 


21.40 


20.86 


11.52 


0.00 


-20.20 




±0.09 


±1.59 


±0.28 




±0.43 


NGC 4593 


20.63 


12.44 


iz.yjL 


U.Ul 


-20.87 




±0.17 


±1.30 


±0.37 




±0.46 


NGC 5548 


21.49 


11.03 


12.89 


0.06 


-21.31 




±0.26 


±3.20 


±0.71 




±0.78 


NGC 7469 


19.66 


5.21 


12.93 


0.06 


-21.15 




±0.14 


±0.44 


±0.32 




±0.45 


Mrk 79 


21.41 


5.68 


14.46 


0.16 


-20.29 




±0.28 


±0.70 


±0.45 




±0.56 


Mrk 279 


20.31 


3.89 


14.53 


0.06 


-20.86 




±0.15 


±0.41 


±0.39 




±0.51 


Mrk 590* 


15.99 


0.71 


14.79 


0.03 


-20.25 




±0.03 


±0.02 


±0.21 




±0.39 


Mrk 817* 


18.25 


0.98 


16.54 


0.44 


-18.58 




±0.09 


±0.13 


±0.39 




±0.51 



^ fie and fio arc the surface brightness fitting parameters for r-'^''^-law 

and exponential forms of the luminosity profile respectively. 

^ re and ro are the scale length fitting parameters for r^/^-law and 

exponential forms of the luminosity profile respectively. 

* indicates an exponential bulge was fitted to the data; all others 

are r^^"^ law fits. 



what might be expected from the uncertainties in the 
measurements . 

4. NOTES ON INDIVIDUAL SEYFERT 1 GALAXIES 

NGC 3227 — The active continuum is weak, and 
the Ca II triplet lines a re seen at nearly full strength 
lINelson fc Whitt]ell200l . Note the presence of Paschen 
absorption lines in figure (21 (particularly P14 A8598.4 be- 
tween the two stronger Ca II triplet lines) indicating the 
presence of a young stellar population. An F5 III star was 
used for a velocity template in this case. The bulge of this 
galaxy is s urprisingly small and faint given its Hubble 
type fSABa ide Vaucouleurs et al]ll99H) a nd , perhaps 
as a result of the dusty nuclear regions iMalkan et al. I 
©98). 

NGC 4051 — Weak Ca II triplet in emission can be 
seen in this object as well as strong and relatively narrow 
01 A 8446 and [Fe II] A 8618 emission. All of these were 
removed by combining model nuclear emission with a 
seeing star as described in section 

NGC 4-151 — Nuclear emission features in the Ca II 
triplet spectrum were removed by combining model nu- 
clear emission with a seeing star. Long and short expo- 
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sure images for this object were obtained with the F547M 
filter. This allowed an accurate subtraction of the nu- 
cleus. Also this filter excludes most of the optical emis- 
sion lines so no correction was applied for the NLR. The 
magnitude in table |21 has been converted to the F606W 
bandpass using CALCPHOT assuming the spectrum is 
well represented by the Sb template in the SYNPHOT 
package. This is the most discrepant point in figure |S1 
for which McLurc & DunloD (2001 ) used the bulge-disk 
decomposition of Baggctt. Bag gett fc AndersonI l)1998fl 
based on ground based photographic imaging. 

NGC 4593 — Nuclear emission features were removed 
by combining model nuclear emission with a seeing star. 

NGC 5548 — Since the velocity dispersion falls off 
quite rapidly with distance from the nucleus a narrow 
spectral extraction was used. 

NGC 7469 — The position angle and longslit extrac- 
tion were c hosen to avoid knot s in the well-known star- 
burst ring llWilson et al.llT99l . 01 A8446 was removed 
with Gaussian fit. 

Mrk 110 — This was the poorest spectrum in our sam- 
ple and as a consequence we quote a large error estimate 
for the velocity dispersion. 01 A8446 removed with Gaus- 
sian fit. 

Mrk 279 — Broad, weak 01 A8446 was removed in 
normal continuum subtraction. 

Mrk 335 — Ca II triplet in emission 
ijMalkan fc FiliDDenkolll98^ could not be removed. 

Mrk 509 — No Ca II triplet features were detected due 
to the strong nuclear continuum. 

Mrk 590 — The bulge of this galaxy is better fit by a 
exponential than an r^^^ law. 

Mrk 817 — 01 A8446 was removed with double Gaus- 
sian fit. The bulge is better fit by a exponential than an 
r^^^ law. 

3C 390.3 — Due to its high redshift, the Ca 11 triplet 
lines fall in a region relatively free of night sky emission 
lines eliminating the need for nod-and-shuffle sky sub- 
traction. However, the atmospheric absorption is severe; 
careful division by atmospheric template was remarkably 
successful. 

5. THE M, - a, RELATION FOR ACTIVE GALAXIES 

In Figure we plot as • symbols our cr* values ver- 
sus Af, values from reverberation mapping. We use the 
M, values directly from the literature, but note that 
these are calculated assuming random orb i ts for BLR 
clou ds following equ ation 5 in IKaspi et all pOOOl) (see 
also lWandel et al. 1 ll999l[) . We discuss the significance of 
this assumption in section|Bl Furthermore, as mentioned 
in the note to table ^ the black hole masses are those 
calculated using the r.m.s. H/3 FWHM rather than the 
mean values, i.e. th e variable component of the line (see 
iWandel et aLlll999|) . Having tried both, we find no sig- 
nificant differences in our results using one or the other. 
For comparison, we also plot as -|- symbol s, elliptical 
and SO galaxies from lGebhardt et al.l l)2000a|) . where Af, 
has been determined from stellar dynamical orbit-based 
modeling. T he solid line is the deriv ation of the Af, — 
relation from lTremaine et alJ l)2002y) and the dashed line 
is a fit to the Seyferts using the o rdinary least squa res bi- 
sector method (hereafter OLSB. ITsobe et aT 1990). The 
fourteen Seyferts, by themselves, show a strong correla- 
tion (linear correlation coefhcient, R = 0.73 for the r.m.s. 




M(F606W) (this paper) 

Fig. 5. — Our Mh,,] measu rements are compared with 
those of IMcLure &: Dunlop I 120011) rescaled to rescaled to Hq = 
75kms-lMpc-l. The soUd hue is X = Y. 



masses with probability of no correlation, Pnuii = 0.29% 
and slightly better for the mean masses R — 0.80, and 
Pnull = 0.065%). 

In figure there appears to be no significant evidence 
for an offset for Seyferts from the iTremaine 'eranil200l 
relation. Computing the residuals of the AGN to the 
solid line we find an average shift of AlogAf,= 0.21 in 
the direction of lower black hole masses at a given cr*, 
with a standard deviation of 0.46 dex for the 14 AGN. 
Clearly, the significance of such a shift is marginal and 
distinctly not in agreement with the factor of 5 found by 
IHoI 1)19991) using the M, - Mb^i relation. 

The exact value of the slope of AT, - cr* has been a 
matter of some debate. The slope, 4.1 ± 0.5 for the AGN 
sample is not significantly different from the slope for 
the normal galaxies, although given the uncertainty, it 
is beyond the ability of this dataset to provide a deci- 
sive result. A larger sample of AGN might well provide 
interesting constraints. 

T hese results ef fec tively co nfirm th e conc lusions 
of IGebhardt et alJ ll2000hD . INelsonI Emh and 
iFerraresTeT^LT l|2001|) that there is no distinct dif- 
ference in the M, — Ue relation between active and 
non-active galaxies. Furthermore, it suggests that 
the reverberation mapping technique for determining 
black hole masses in AGN produces results which are 
consistent with masses determined by stellar and gas 
dynamical techniques. We should emphasize, however, 
that this is an indirect check on reverberation mapping. 
A direct verification may be possible should a target 
of reverberation mapping campaigns go into a low 
luminosity phase, allowing application of standard 
stellar dynamical analysis. Indeed an HST cycle 12 
proposal involving target of opportunity observations 
of NGC 3227 and NGC 4151 has been approved (Prop. 
No. 9849, PI B. Peterson). 
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Fig. 6. — Black hole mass is plotted against velocity dispersion for 
Seyfert galaxies from Table 111 shown as » symbols an d for normal 
elliptical and SO galaxies from lGebhardt et al. f2000a'). shown as + 
symbols. The black hole masses for Seyfert galaxies arc determined 
from reverberation mapping, while the tjebhardt ct al. |20003) 
masses are from dynamical modeling of the nuclear stellar kine- 
matics. The solid line is the relation calculated bv iTremaine et al.l 
J2002) and the dashed line is a fit to the Seyferts using the ordinary 
least-squares bisector method llsobe et aLlllflfllTl . 



Our results, of course, are very encouraging for inves- 
tigations of the role of black hole mass in AGN. How- 
ever, the sample used here is by necessity composed of 
lower luminosity Seyfert 1 galaxies. Further study of the 
M, — (Te relation is needed to validate its application to 
more luminous AGN. 

6. THE M. - Mbui RELATION FOR AGN REVISITED AND 
THE FABER-JACKSON RELATION 

It is essentially conventional wisdom that active nuclei 
are preferentially found in bulge dominated galaxies (e.g. 
Wfi- Filippcnko & Sargent 1997). Confirmation of the 
M, — cr» relation for Seyferts is important support for this 
tenet of AGN physics. It strongly reinforces the notion 
that the bulge gravitational potential plays an important 
role in determining the properties of activity in galaxies. 
However, if we accept that the M, — cr, relation is the 
same for Seyferts as it is for normal bulges and ellipticals, 
what do we make of previous claims of smaller A/, values 
in AGN for a given bulge luminosity, discussed in the 
introduction to this paper? 

In Figure [3 we plot Mi,ui for the Seyfert 1 galaxies 
in table 121 versus the reverberation mapping M, values. 
Again for comparison, ellipticals and bulges from the 
iGebhardt et al.l l)2000a(l sample, are plotted as -I- sym- 
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Fig. 7. — Black hole mass is plotted against bulge absolute 
magnitude for Seyfert galaxies (• symbols) and the Gebhardt et a'Q 
iSOOOa) sample (+ symbols) rescaled to Hq = 75kms^-'^Mpc^^. 
The solid line is a fit to the elliptical galaxies. 



bols, where we have converted their abs olute V magni- 
tudes as pubhshed in iFaber et ail l|1997D to the F606W 
bandpass and scaled to Hq = 75 kms~"'^Mpc~ ''^. Given 
the narrow range of colors for the lGebhardt et a l. ( 20004 
sample {{B ~ V) ^ 0'?95, std. dev. 0™047) and the 
expectation that the Seyfert bulge colors are likely to 
be more scattered, we choose to apply the correction to 
the ellipticals. The conversion was determined using the 
template elliptical galaxy from the SYNPHOT package 
and the CALCPHOT task (Mpeoew = V"-0'"33). 

Although not a strong trend, the majority of Seyferts 
lie below and to the right of the OLSB fit to the ellipti- 
cals, sho wn as the s olid line. This is consistent with the 
report bv lHol l)1999D of lower black hole masses at a given 
bulge magnitude. The median deviation in black hole 
mass for the Seyferts from the fit to the normal galax- 
ies is 0.4 dex, or about a factor of 2.5. There is really 
no significant correlation for the Seyferts, (linear corre- 
lation coefficient R = —0.16), compared with a strong 
trend found for the ellipticals {R = —0.80). This may 
be due in part to the smaller sample size, but this can- 
not be the whole story since the Seyferts show a strong 
correlation between M, and cr* (see section . 

We suggest that it is the bulge luminosities, not the 
black hole masses, which are the origin of the deviation. 
Taken in this way the bulges of Seyferts are on average 
0™7 brighter than normal ellipticals with similar black 
hole masses based on the distribution of points in figure 
Strong support for this interpretation is found by 
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1.8 2 2.2 2.4 2.6 

Log (7 (km/s) 

Fig. 8. — The Faber-Jackson plot for Seyfert 1 galaxies 
from Table p| symbo ls) and for elliptical and SO galaxies from 
IGebhardt etall 12000al. + symbols). The solid line is a fit to the 
elliptical sample and the dashed line is a fit to the Seyfert galaxies. 
The Seyfert galaxies in the sample show a tendency for brighter 
Mjuj at a given cr» . 



plotting the Faber-Jackson diagram, tr* vs. Mbui, for 
both samples, shown in fig. |S| We find that that at a 
given cr* the Seyfert galaxies tend to lie above the fit 
to the normal galaxies. The mean residual in Mb^i for 
the Seyferts relative to the normal galaxy fit is —0^4. 
Furthermore there is a moderately strong correlation for 
the Seyferts, {R = —0.6). This improvement over the 
results for Mi,ui vs. M, naturally suggests that Mt,ui 
is more closely related to ct* than to Af,, and that the 
increase in bulge luminosities with black hole masses is 
secondary and follows through their mutual trends with 
velocity dispersion. 

The result is identical to that of iNelson fc Whittf^ 
(^90) who found a similar offset for a much larger sam- 
ple, consisting of mostly Seyfert 2s . The r e sult al so 
confirms the earlier conclusions of iWhittld l)1992bl Ol 
who used galaxy rotation amplitude and FWHM[oiii] as 
virial parameters, to demonstrate a tendency for brighter 
bulges in Seyferts than normal galaxies of the same bulge 
mass. Some ambiguity remained, however, since the 
bulge magnitudes in those studies were not obtained by 
bulge/disk decomposition. Instead, after correcting for 
nuclear emission, the relation between bulgc-to-disk ratio 
and Hubble stage of Simicn & dc Vaucoulcurs (1986) was 
used to estimate bulge magnitudes. This t echn ique has 
been c riticized bv McLure fc Dunlop I l)200lD and lWandel 
1)2002(1 who suggest that it results in overestimates of 



Mbui ■ The present confirmation of the offset in the Faber- 
Jackson relation for Seyferts using explicit bulge-disk- 
nucleus separation, however, suggests that while there is 
admittedly significant scatter, any systematic errors are 
small. 

In contras t to thelHnl lll99 9D result. iMcLure fc Dunlop I 
||2001D and iWandell l| 2002l claim no deviation in the 
M, — Mbui relation for Seyfert and normal galax- 
ies. We can reconcile th is conflict by recognizing that 
iMcLure fc Dunlopi f200l) choose to calculate M, assum- 
ing a disk-like morphology and purely rotational kine- 
ma tics for the BLR as opposed to random orbits assumed 
bv iKaspj et alJ l|200(]D . Thus the observed line widths 
depend on the orbital velocity, V, and the inclination 
of the BLR disk. After making assumptions about the 
range and mean of the distribution of inclinations, they 
adopt a relation between orbital velocity and the BLR 
line widths, V = l.bx FWHM, distinctly different from 
the case for random cloud orbits, V = x FWHM. 

This results in M, values that are a factor of 3 larger 
than the standard reverberation mapping calculations in 
which the BLR velocity field is assumed to be random. 
Therefore, if we were to couple their bulge measurements 
with the M, values calculated assuming random BLR 
motions, we should also get brighter bulges at a given 
M, for AGN. Of course the specific form of the BLR ve- 
locity field is unknown. Our results, at least for the case 
of random cloud motions, are consistent with an identi- 
cal A/, — (Je relation for Seyfert and normal galaxies, and 
a tendency for brighter bulges in Seyfert s. 

A different argument applies to the iWandel l|2fl02l) 
analysis, which assumes random motions dominate in 
the BLR. This study reports a correlation between H/3 
line width and the suspected o verestimate of Mh,d result- 
ing from the application of the lSimien fc de VaucouleursI 
^986) B/T -morphology relation discussed above. This 
correlation is then used to correct the Mh^i values, re- 
sulting in a reduction of the bulge luminosities for much 
of the sample. However, the physical link between BLR 
kinematics and differences in B/T is not established. It 
is therefore unclear why such a correction should be re- 
quired. Removing this correction from Wandel's results 
should also result in relatively brighter bulges for AGN. 

The most straightforward interpretation of these re- 
sults is that Seyfert galaxies have, on average, brighter 
bulges as a result of lower M/L ratios. Thus, the differ- 
ence is one of stellar population, indicating that Seyferts 
are more likely than non-active galaxies to have had re- 
cent episodes of nuclear star formation. This is certainly 
plausible, given the voluminous research on the link be- 
tween AGN and starburst galaxies. It is informative to 
use model stellar populations to estimate the age and 
mass fraction of the bulge required to produce a brighten- 
ing of this degree . Guy Worthey's Dial-a-Galaxy website 
ijWorthev 1 11994') allows users to mix stellar population 
models of different ages, metallicities, etc., to estimate 
the changes to the integrated light properties of galaxy 
spectra. A simple run, with two solar metallicity popula- 
tions of 1 Gyr and 12 Gyr, shows that if 15% of the total 
mass is from the young population the i?-band magni- 
tude will be ^ 0™5 brighter than a population composed 
of completely old stars. Obviously more complex and 
realistic models could produce similar results. Neverthe- 
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less, this simple estimate indicates a typical size and age 
for a burst of star formation to produce the shift seen in 
the Seyfert sample. 

It is interesting to note that a number of recent 
morphological studies of Seyfert galaxies have found 
no evidence for morphological triggers of nuclear star- 
bursts — galaxy inte ractions and bars, including sub 



kpc scale bars (e.g. 'Ho. Fihppenko fc Sargent I Il997t 
IPogge fc Martiml 12002: Regan & Mulchaev 2001). Al- 
though it is conceivable that there is a time delay be- 
tween a galaxy interaction and the onset of nuclear fu- 
eling, the result is troubling since we might expect at 
least some evidence of the dist urbance to remain. A pos- 
sible solution was hinted at in lNelson fc Whittld l|T99l 
based on the fact that Seyferts are not found outside 
the general envelope of scatter for normal galaxies in 
the Faber- Jackson plot. Considering the residuals to the 
Faber-Jackson relation, they found that the distribution 
of A log (J* is considerably narrower yet overlapping that 
for non-active galaxies, a detail which is reproduced in 
Figure |H1 Assuming that the spheriods of normal galax- 
ies exhibit a broad range of star-formation histories, per- 
haps Seyferts are more akin to normal spirals with above 
average central star formation rates than to massive star- 
bursts initiated by strong interactions or bars. 



7. SUMMARY 

We have obtained new measurements of the bulge stel- 
lar velocity dispersion in Seyfert 1 galaxies with nuclear 
black hole masses determined using the reverberation 
mapping technique. Our results have shown that Seyfert 
galaxies follow the same correlation between M, and cr, 
as for non-active galaxies and that reverberation map- 
ping yields black hole masses which are consistent with 
masses determined through stellar dynamical modeling 
techniques. We also investigated the origin of previous 
assertions that the masses of black holes in active galaxies 
were smaller than those in normal galaxies (or at least 
underestimated) at a given bulge luminosity. We find 
that the best explanation is that the bulges are brighter 
on average in Seyfert galaxies, suggesting lower M/L ra- 
tios, most likely as a result of higher star formation rates 
in their recent history. 
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